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Abstract 

We present a model for leptonic mixing and the lepton masses based on flavor 
symmetries and higher-dimensional mass operators. The model predicts bilarge lep- 
tonic mixing (i.e., the mixing angles 612 and #23 are large and the mixing angle 
$13 is small) and an inverted hierarchical neutrino mass spectrum. Furthermore, it 
approximately yields the experimental hierarchical mass spectrum of the charged 
leptons. The obtained values for the leptonic mixing parameters and the neutrino 
mass squared differences are all in agreement with atmospheric neutrino data, the 
Mikheyev-Smirnov-Wolfenstein large mixing angle solution of the solar neutrino 
problem, and consistent with the upper bound on the reactor mixing angle. Thus, 
we have a large, but not close to maximal, solar mixing angle 612, a nearly maximal 
atmospheric mixing angle 623, and a small reactor mixing angle #13. In addition, 
the model predicts 612 — \ — #13- 
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1 Introduction 



The fermionic sector of the standard model (SM) of elementary particle physics 
is described by 13 renormalized parameters (6 quark masses, 3 charged lep- 
ton masses, 3 CKM mixing angles^, and one CP violation phase). Obviously, 
these parameters are not arbitrary, but exhibit relations which can only be un- 
derstood when going beyond the SM. One possibility to obtain realistic quark 
masses and CKM mixing angles in an extension of the SM is to introduce 
flavor symmetries that are sequentially broken. At the first glance, the hier- 
archical mass spectra of the quarks and the charged leptons actually suggest 
underlying non-Abelian flavor symmetry groups acting on the first and second 
generations only.0 However, in the light of recent atmospheric [5-9] and so- 
lar [10-17] neutrino experimental results, it seems to be difficult to extend this 
idea to the neutrinos. Especially, the result that, among the different possible 
solutions of the solar neutrino problem, the Mikheyev-Smirnov-Wolfenstein 
(MSW) [18-20] large mixing angle (LMA) solution is the presently preferred 
one [21-23]Q, draws a picture of the involved flavor symmetries and their 
breaking mechanisms that differs remarkably from the early approaches, which 
have been applied to the quark sector. In the "standard" parameterization, 
the MSW LMA solution implies that we have a bilarge mixing scenario in the 
lepton sector in which the solar mixing angle 612 is large, but not necessarily 
close to maximal, the atmospheric mixing angle #23 is nearly maximal, and 
the reactor mixing angle # 13 is small. Clearly, this is in sharp contrast to the 
quark sector in which all mixing angles are small [28] and it indicates that the 
flavor symmetries act on the third generation too. 

By assuming only an Abelian U(l) flavor symmetry, one obtains that the at- 
mospheric mixing angle may be large, but cannot be enforced to be nearly 
maximal [29-31]. Therefore, a natural close to maximal z/ A1 -z/ r -mixing can be 
interpreted as a strong hint for some underlying non-Abelian flavor symme- 
try acting on the second and third generations [32,33]. Neutrino mass models 
which give large or maximal solar and atmospheric mixing angles by putting 
the second and third generations of the leptons into the regular representa- 
tion of the symmetric group S2 [34] or into the irreducible two-dimensional 
representation of the symmetric group S3 [35] are, in general, plagued with 

3 The mixing angles in the quark sector are usually called the Cabibbo-Kobayashi- 
Maskawa (CKM) mixing angles [1,2]. 

4 By placing the first two generations into irreducible representations of flavor sym- 
metries, one can in supersymmetric models achieve near degeneracy of the corre- 
sponding squark masses thus suppressing large flavor changing neutral currents [3,4]. 

5 Actually, there exist several recent global solar neutrino oscillation analyses in- 
cluding the latest SNO data that strongly favor the MSW LMA solution of the solar 
neutrino problem, see, e.g., Refs. [16,24,25,23,26,27]. However, we have chosen to 
list and use the values obtained in Ref. [23]. 
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a fine-tuning problem in the charged lepton sector, since they tend to pre- 
dict the muon and tau masses to be of the same order of magnitude, i.e., 
they lack providing an understanding of the hierarchical mass spectrum in the 
charged lepton sector. A recently proposed model based on an SU(3) flavor 
symmetry [4] gives approximately bimaximal leptonic mixing as well as a suc- 
cessful description of the charged fermion masses, but predicts the presently 
disfavored MSW low mass (LOW) or vacuum oscillation (VAC) solution of 
the solar neutrino problem. Similarly, the highly predictive models of flavor 
democracy [36-43] yield large solar and atmospheric mixing angles, but fit 
the LOW or VAC solution rather than the MSW LMA solution [44]. In grand 
unified theory model building, it seems that the MSW LMA solution with 
a normal hierarchical neutrino mass spectrum is more natural than with an 
inverted one [45] (for a recent phenomenological analysis of minimal schemes 
for the MSW LMA solution with inverted hierarchical neutrino mass spectra, 
see, e.g., Ref. [46]). Also from empirical lepton and quark mass spectra anal- 
yses a normal hierarchical (or inverse hierarchical) neutrino mass spectrum 
seems to be rather plausible [47] . A comparably simple way of generating the 
MSW LMA solution with normal hierarchical neutrino mass spectra is, e.g., 
provided by models based on single right-handed neutrino dominance [48]. 

In a previous Letter [49], we introduced a model for bilarge leptonic mixing 
based on higher-dimensional operators, using the Froggatt-Nielsen mecha- 
nism, and Abelian horizontal flavor symmetries of continuous and discrete 
types. In this paper, we consider a modified and extended version of this 
model and we explicitly demonstrate the vacuum alignment mechanism, which 
produces a nearly maximal atmospheric mixing angle 823 as well as a large, 
but not close to maximal, solar mixing angle 612, as required by the MSW 
LMA solution. Simultaneously, the vacuum alignment mechanism generates 
a strictly hierarchical charged lepton mass spectrum thus resolving the fine- 
tuning problem many realistic models, which seek to predict the MSW LMA 
solution, are suffering from. Furthermore, this model gives a small mixing of 
the first and second generations of the charged leptons, which is comparable 
with the mixing of the quarks, whereas the mixing among the neutrinos is 
essentially bimaximal [50]. The actual leptonic mixing angles are then a result 
of combining the contributions coming from both the charged leptons and the 
neutrinos (see Fig. 1). Thus, the model predicts the relation 9 12 — | — #13 
between the solar mixing angle 612 and the reactor mixing angle #13, which is 
non-zero and lies in the range of the quark mixing angles. 

Note that our study assumes that there are three neutrino flavors, and there- 
fore, three neutrino flavor states u a (a = e, \i, r) and also three neutrino mass 
eigenstates v a [a = 1,2,3). Furthermore, it assumes that all CV violation 
phases are equal to zero. 

The paper is organized as follows: In Sec. 2, we introduce the representation 
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Fig. 1. Flow-scheme for computing leptonic mixing parameters and lepton masses 
from any given charged lepton and neutrino mass matrices. 

content of our model including U(l) charges and discrete symmetries. The 
multi-scalar potential of the model is then analyzed in Sec. 3, where we ex- 
plicitly demonstrate the vacuum alignment mechanism. Next, in Sees. 4 and 5, 
the Yukawa interactions of the charged leptons and neutrinos, respectively, are 
investigated and discussed, which lead to the mass matrices of the correspond- 
ing particles. In Sec. 6, the lepton mass matrices are diagonalized yielding the 
charged lepton masses and the neutrino mass squared differences as well as 
the charged lepton and neutrino mixing angles (see again Fig. 1). In Sec. 7, 
the total leptonic mixing angles are derived and calculated. Implications for 
neutrinoless double /3-decay, astrophysics, and cosmology are briefly studied 
in Sec. 8. Finally, in Sec. 9, we present a summary as well as our conclusions. 
In addition, we present in the Appendix a scheme for how to transform any 
given 3x3 unitary matrix to the "standard" parameterization form of the 
Particle Data Group [28]. 



2 The representation content 



Let us consider an extension of the SM in which the lepton masses are gener- 
ated by higher-dimensional operators [51,52] via the Froggatt-Nielsen mecha- 
nism [53]. (A classification of effective neutrino mass operators has been given 
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in Ref. [54].) Since we are mainly concerned with the question of whether there 
is a possible naturally maximal z/ M -z/ T -mixing in the MSW LMA solution or 
not, for which the properties of the quarks are seemingly irrelevant, we will, for 
simplicity and without loss of generality, omit the quark sector in our further 
discussion. For a recent related study including also the quark sector, see, e.g., 
Refs. [55,56]. In a self-explanatory notation, we will denote the lepton doublets 

as L a = {v aL e a i)i where a = e,\i,r, and the right-handed charged leptons 
as E a = e a R, where a = e,fi,r. The part of the scalar sector, which carries 
non-zero SM quantum numbers, consists of two Higgs doublets Hi and H 2 , 
where Hi couples to the neutrinos and H2 to the charged leptons. This can be 
achieved by assuming, e.g., a discrete Z 2 symmetry under which H 2 and E a 
(a = e, /i, r) are odd and Hi and the rest of the SM fields are even. The masses 
of the charged leptons arise through the mixing with additional heavy right- 
handed charged fermions, which all have masses of the order of magnitude of 
some characteristic mass scale Mi. Apart from some general prescriptions of 
their transformations under the flavor symmetries, which will be introduced 
below, it is not necessary to explicitly present the fundamental theory of these 
additional (or extra) charged fermions. This is in contrast to the neutrino sec- 
tor, which we will extend by five additional heavy SM singlet Dirac neutrinos 
N e , Np, N T , Fi, and F 2 . The neutrinos Fi and F 2 are supposed to have masses 
of the same order Mi as the charged intermediate Froggatt-Nielsen states, 
whereas N e , N^, and N T all have masses of the order of magnitude of some 
relevant high (unification) mass scale M 2 . While M 2 takes the role of some 
seesaw scale [57-59] (and it is therefore responsible for the smallness of the 
neutrino masses), Mi can be as low as several TeV [60,61]. In order to ob- 
tain the structures of the lepton mass matrices from an underlying symmetry 
principle in the context of a renormalizable field theory, we will furthermore 
extend the scalar sector by additional SM singlet scalar fields <pi, <p 2 , • • • , 0io, 
(p'i,(p' 2 , ■ ■ • ,0' 6 , and 9 and we will assign these fields gauged horizontal U(l) 
charges Qi, Q 2 , and Q3 as follows: 



(Q1.Q2.Q3) 



iV M .JV x 
Fi 




(1,0,0) 
(0,1,0) 
(1,0,0) 
(0,1,0) 
(1,0,0) 

(-1,0,1) 



Hi 

4>i 



03 
05 

<t>' 3 

07 



09 
010 



e 



04 
06 




(0,0,0) 
(1.-1.2) 
(0, 0, 0) 
(0, 0, 1) 
(0, 0, 0) 
(-1,-1,0) 

(-2,0,1) 
(0, 0, 0) 

(0,0,-1) 
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In the rest of the paper, it will always be understood that the Higgs dou- 
blets Hi and H 2 are total singlet under transformations of other additional 
symmetries. Note that the charges Qi and Q 2 are anomalous. However, it is 
known that anomalous U(l) charges may arise in effective field theories from 
strings. Then, the cancellation of the anomalies must be accomplished by the 
Green-Schwarz mechanism [62]. 

The first generation of the charged leptons is distinguished from the second 
and third generations if we require for P = e 27rl / n , where the integer n obeys 
n > 5, invariance of the Lagrangian under transformation of the following Z n 
symmetry: 



E P . 



Vi:l 



P- 4 E £ , -> E T -> p-'Er, 



01 -> P(f>l, 02 -> -P02, 
03 -> P(f>3, 04 -> -P04, 
05 -> ^05, 06 -> -P06, 



where we assume that the fields iV e , A*^, iV T , Fi, and F 2 are singlets under 
transformation of the symmetry T>i. In addition, the symmetry T>i forbids 
the fields 0i, 02, . . . , 06 to participate in the leading order mass terms for the 
neutrinos. Furthermore, the permutation symmetries 



0i <"»• 02, 01 ^ 02, 
07 -> -07, 



(2) 



X>4 : < 



03 ^ 04, 03 <"»• 04, 
05 ^ 06, 05 «"»• 06, 
07 ~> -07, 

02 -> -02, 04 -> -04, 06 
07 08 



(3) 



-06, 



(4) 
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are responsible for generating a naturally maximal atmospheric mixing angle, 
since they establish exact degeneracies of the Yukawa couplings in the leptonic 
2-3-subsector. These permutation symmetries also play a crucial role in the 
scalar sector in which they restrict some of the couplings in the multi-scalar 
potential to be exactly degenerate (at tree level), which means that degenerate 
vacuum expectation values (VEVs) can emerge after spontaneous symmetry 
breaking (SSB). This so-called vacuum alignment mechanism can work if we 
assume the discrete symmetries 



0'i - -<; 



3> ^3 
5» 05 



-01, 
-03, 
-05, 



P- k 4>'i, 02 



03- 
05" 
01 
03 

05 

09 - 



* P fe 01, 02 



> P'03, 
■ P™05, 



^4 ■ 

iio 



(5) 



PiV P 



P 



y 2, 



' ^04, 
P~V 6 , 
* P fc 02, 
- P'04, 
► P m 06, 
"> P01O, 



(6) 



where /c, /, and m are some integers. For the symmetry V 6 we additionally 
require that the Froggatt-Nielsen states with non-zero hypercharge can only 
be multiplied by factors P™, where n is an integer multiple of k, I, or m, 
and that the differences \k — 1\, \k — m\, and \l — m\ are sufficiently large. The 
only fermion with non-vanishing hypercharge that transforms differently is the 
right-handed electron E e . These symmetries restrict the allowed combinations 
of the scalar fields in the higher- dimensional lepton mass operators as well 
as in the renormalizable terms of the multi-scalar potential. Thus, possibly 
dangerous terms in the multi-scalar potential can be forbidden, which could 
otherwise spoil the vacuum alignment mechanism. 



At this stage, it is appropriate to point out some implications concerning the 
nature of the discrete symmetries. It has been found that relations between 
Yukawa couplings established by standard discrete symmetries can only re- 
main unbroken by quantum gravity corrections if the discrete symmetries are 
gauged [63]. These "discrete gauge symmetries" appear in continuum theo- 
ries when a gauge symmetry group G is broken to a discrete symmetry sub- 
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group H. Since acceptable continuous gauge theories have to be free from 
chiral anomalies, one obtains discrete anomaly cancellation conditions, which 
strongly constrain the massless fermion content of the theory [64]. At a more 
fundamental level, this implies for our model that the permutation symme- 
tries must actually be gauged, since the relations, which are based on them, 
will prove to be crucial for obtaining essentially strict maximal atmospheric 
mixing. However, it is interesting to note that in M-theory, one expects the 
discrete symmetries to be always anomaly-free [65]. 



3 The multi-scalar potential 

3. 1 The two-Higgs doublet potential 

From the most general two-Higgs doublet potential of the fields H\ and H 2 
(for an extensive review on electroweak Higgs potentials, see, e.g., Ref. [66]) we 
conclude that, in presence of the Z 2 symmetry, which distinguishes between Hi 
and H 2 (see Sec. 2), each of these fields appears only to the second or fourth 
power in the multi-scalar potential. Hence, in any renormalizable terms of 
the multi-scalar potential, which mix Hi or H2 with the SM singlet scalar 
fields, the Higgs doublets are only allowed to appear in terms of their absolute 
squares \Hi\ 2 and \H 2 \ 2 . Next, since the Higgs doublets carry zero Qi,Q 2 , 
and Q 3 charges and are Dj-singlets, where i = 1,2,... ,6, there exists a range 
of parameters in the multi-scalar potential for which the standard two-Higgs 
electroweak symmetry breaking is possible. Furthermore, this implies that we 
can, without loss of generality, separate the SM singlet scalar part from the 
Higgs-doublet part in the multi-scalar potential by formally absorbing the 
absolute squares of the VEVs \{Hi)\ 2 and \ (H 2 }\ 2 into the coupling constants 
of the mixed terms. Then, since the vacuum alignment mechanism of the SM 
singlet fields is independent from the details of the Higgs doublet physics, 
we can in what follows discard the effects of the Higgs doublets and fully 
concentrate on the properties of the SM singlet scalar fields. 

3.2 Interactions of the fields 9 and (f> w 

The symmetry T> 6 requires that the fields 9 and (f) W enter the renormalizable 
interactions of the scalar fields only in terms of the operators 0g0io, 09</>Io> |</>9| 2 > 
and |</>io| 2 - The product 0g0io has the U(l) charge structure (—2,0, 1), imply- 
ing that the only renormalizable interaction in the scalar potential, which 
involves this product, is actually ~ |0 9 | 2 |0i O | 2 . We can assume that both of 
the fields 0g and (fi w finally develop non-vanishing VEVs with magnitudes 
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|(09)| and |(0io)| that are not too large. Since the fields 0g and 0io are singlets 
under transformations of all the permutation symmetries T> 2 , T> 3 , and P4, they 
will have no effect on the relative alignment of the rest of the scalar fields, be- 
cause they enter the corresponding interactions only in terms of the absolute 
squares |(0g)| 2 and |(0i O )| 2 . For this reason we will, without loss of generality, 
discard the terms in the scalar potential which involve the fields 0g and 0io in 
our considerations. 

From the assignment of the U(l) charges Qi,Q 2 , and Q 3 and the symme- 
try T>§ (which does not permute any of the fields) it follows that any renor- 
malizable term in the scalar potential which involves the SM singlet fields 
01)02,05,06,07,08, or 9 can only be allowed if these fields appear in one of 
the following combinations: 

01,201,2 , 05,605,6 , 05,6# , 07,807,8 » > ( 7 ) 

where "0ij" (i, j = 1, 2, . . . , 8) denotes some linear combination of the fields 0^ 
and (pj. Among the products in Eq. (7) only the product 05,6# transforms non- 
trivially under the symmetry T>\. In addition, since the fields 0^,02,... , 0' 6 
are Pi-singlets, we observe that the fields 03 and 04, which carry a non- 
zero Di-charge, can only appear either in the combination 03,403 4 or in the 
combination 03 j4 05,6#- However, the latter combination is forbidden by XV 
invariance. 

3. 3 Interactions of the field 9 

As in Sec. 3.2, we can assume that the field 9 finally develops a non-vanishing 
VEV with a magnitude \(9)\ that is not too large. Except for the combination 
05,6$ in Eq. (7) all scalar interactions involve an equal number (0, 1, or 2) of 
the fields 9 and its adjoint 9\ which can then be paired to the absolute square 
|#| 2 . Since the interaction 05,6#|#| 2 is forbidden by the symmetries V\ and 
Vq and 9 is a total singlet under transformations of the discrete symmetries 
T>i,T>2, ■ ■ ■ ,T> 6 , all terms which involve the absolute square \9\ 2 will have no 
influence on the relative alignment of the rest of the scalar fields. For this 
reason we can, without loss of generality, omit the terms in the scalar potential 
which involve \9\ 2 in our considerations. 

3.4 Interactions of the fields 7 and 8 

From the U(l) charge assignment we conclude that only an even number of 
the fields 07 and 0s (or their complex conjugates) can participate in the scalar 
interactions. Let us now especially consider the operator 7 0g (or equivalently 
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its complex conjugate). Under application of each of the symmetries T>2 and 
£>3 the operator 7 08 changes sign. Therefore, the symmetry T> 2 requires this 
operator to couple to one of the fields taken from the set {0^, <f>' 2 , 0i, 2 , 07}- Si- 
multaneously, the symmetry V 3 requires this operator to couple only to fields 
taken from the extended set {0' 3 , 0' 4 , 0' 5 , 0' 6 , 03, 04, 05, 06, 07}. Both conditions 
can only be fulfilled if the operator 0708 couples to the field 07 or its com- 
plex conjugate, but not to the operator products (07) 2 , (0 7 ) 2 , or |07| 2 . Since 
the product 7 08 carries the U(l) charges (0, 0, 0), the operator 7 0s can only 
couple to some linear combination of the operators 070§ and 0807- As a conse- 
quence, if a general operator 07,807,8 enters an interaction with scalars, which 
are different from the fields 07 and 0s, then this operator 7 , 8 7)8 is a linear 
combination of the absolute squares |0 7 | 2 and |0s| 2 - 

Let us denote by 0j and 4>j two scalar fields, which are different from the fields 
07 and 8 . Taking Eq. (7) and the operator 03,403, 4 into account, the operator 
product 0j0j must be of one of the following types: 

01,201,2 , 03,403,4 , 05,605,6 , 01,20'l,2 > 03,4^3,4 > 05,6^5,6 > ( 8 ) 

where the last three combinations follow from the IVinvariance. Next, the 
symmetry T>^ implies that <pj = 0f for the combinations in Eq. (8), i.e., the 
product is the absolute square 0j0j = |0;| 2 . Using the result of the previous 
paragraph, invariance under transformation of the symmetry D 4 gives for the 
most general interactions of the fields 07 and 8 with the other scalar fields 
the terms 

(I07| 2 + |0s| 2 ) E C *H 2 , (9) 

where (fi can be any of the scalar fields, which are not identical with the fields 
07 or 08 and q are some real-valued coupling constants. (Dimension-three 
terms of the types \<p7\ 2 <pi or |08| 2 y?i, where y?j ^ 07, 08, are forbidden by the 
U(l) charge assignment and the symmetry D 6 , which does not permute any 
fields.) Taking everything into account, the U(l) charge assignment and the 
symmetry £> 4 restrict the most general terms in the scalar potential, involving 
the fields 07 and 0s, to be 



U 7l 8=/i 2 (i07r+i08r)+^(i07i 2 +i08i 2 ) 2 +(i07r+i08r) e c *n 2 

<Pi¥=4>7,<p8 

+ &|07| 2 |0 8 | 2 , (10) 

where /i 2 , re, a, and b are real- valued constants. Parameterizing the VEVs of 



+ a 



(070s) 2 + 
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07 and 0s as 



' (07) ^ 




( e 1/?1 cos a ^ 


- 




v (0s) J 




y e 1/32 sin a j 



where i> is some real- valued number, a is the angle which rotates (07) and (0s), 
and f3\ and /?2 denote the phases of the VEVs, we observe in Eq. (10) that 
the first three terms exhibit an accidental U(l) a symmetry. This symmetry is 
broken by the last two terms, which are therefore responsible for the alignment 
of the VEVs. Note that the scalar potential is symmetric under the exchange 
07 «— > 0s and that the term with coefficient a is ~ K[(070 8 ) 2 1. In Eq. (10), we 
will choose k > and assume the rest of the coupling constants to be negative. 
Then, the lowest energy state is characterized by a = | and (5\ — fii G {0, tt} 
or equivalently 



(0 7 ) = ±(0 8 ), 



(12) 



i.e., the VEVs are degenerate up to a sign. When considering the Yukawa in- 
teractions of the neutrinos, it will turn out that the degeneracy of the VEVs is 
responsible for a nearly maximal atmospheric mixing angle. Thus, we can from 
now on restrict our discussion of the scalar potential to the fields 0i, 02, . . . ,06, 
and <f>[, 2 , ... , 0q. This discussion will follow in the three coming subsections. 



3.5 The potential of the fields 0i, 2 , . . . ,06 



In all two- fold and four-fold products involving only the fields 0j (i — 1, 2, . . . , 6), 
the discrete symmetry V 5 requires the number of these fields, which are de- 
noted by even (or odd) indices, to be even. In the scalar potential, linear and 
tri-linear terms of the fields 0i, 02, . . . ,06 are forbidden by the symmetry T>q. 
Taking the combinations in Eq. (7) and 03,403 4 into account, the allowed two- 
fold products of the fields 0i, 2 , . . . ,06 can only be of the type |0j| 2 , i.e., they 
must be absolute squares of the fields. Similarly, we obtain that all four-fold 
products of the fields 0i, 02, . . . ,06 must be of the types 



(010^) 2 , 01020304 , 01020403 , 01020506 , 01020605 > 

(0 3 0l) 2 , 03010506, 03010605, (050 6 ) 2 , \<t>i?\<t>j\\ \<t>i\\ (13) 

and their complex conjugates, where i,j = 1,2,... ,6. A general four-fold 
product of the types in Eq. (13) can be written as 
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{o4>\4>j + b(f>)(t>i)(t>\(t>i + {c<t>\(t)j + d(f>)(t>i)(t>\(t>k + h.c. 
= (a + d*)4>\4>3 + (b + c*)4>)4>] 4<Pi + [(c + b*)4>\4>j + (d + a*)$&] 0*0 fc , 

(14) 

where a,b,c, and d are complex- valued constants. Assume that i ^ j. Then, 
from Eq. (13) we observe that k ^ I and we can, without loss of generality, 
assume that the index-pairs and (k,l), respectively, combine the fields 
which are interchanged by the discrete symmetry T> 2 or T> 3 . (If i = j, then 
it follows from Eq. (13) that k — I, which will be discussed below.) Let, in 
addition, {i,j} ^ {k, I}. Then, application of the symmetries V 2 and V 3 yields 
a + d* = d + a* and b + c* = c + b*. We can therefore rename the constants as 
a + d* — > a and b + c* — > b, where now a and b are real constants, and write 
the term in Eq. (14) as 



(a0l0j + b4>)(pi)4>l(j)i + (tylfa + a4>)<pi)4>\<p k 

^■(fyj + (t>)(t>i){(t>i(t>i + 4>](f>k) + °^{(t>\(t>3 - (t>)(t>i){(t>{(t>i - <f>\(i>k) 



a + b 



a — b 



3(0L>,). 



(15) 



Since the fields 3 , 4 , 5 , and 06 are singlets under transformation of the 
discrete symmetry T> 2 , we can have a ^ b in the case that ((f>i,<f>j) = (03,04) 
and (0^,0;) = (05, 6 ). However, if (i,j) = (1,2), then application of the 
discrete symmetry T> 2 further constrains the constants in the above general 
form to fulfill a = b, and therefore, the last term in Eq. (15) vanishes. 

As a cause of the symmetries V 2 and V 3 , the products {4>i4>]) 2 in Eq. (13), 
where = (1,2), (3,4), (5,6), appear in the potential always as 



a 



(0,0|) 2 + {4>W 



3 ' ' 

where a is some real- valued constant. 



2a3? 



(16) 



Let us now turn the discussion to the terms |0;| 2 |0fc| 2 in Eq. (13), where % ^ k. 
Assume that the fields 0j and 0& cannot be combined into one of the pairs 
(0i, 02), (03, 04), or (05, 6 ). Then, a general term of this type is on the form 

(a|0,| 2 + &|0,| 2 )|0 fc | 2 + (c|0,| 2 + d|0/)|0*| 2 , (17) 

where a, b, c, and d are real-valued constants and 

(0i, 0j), (0fc, 0/) G {(1, 2), (3, 4), (5, 6)}. Application of the symmetries V 2 and 
£>3 yields the conditions a — d and b = c, and thus, we can rewrite the above 
part of the potential as 

a -^m 2 + mm 2 + n 2 ) + ^m 2 - \mm 2 - i0d 2 )- (is) 
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If = (3,4) and (k,l) = (5,6), then in general a ^ b, since the fields 

03,04,05, and 6 are X> 2 -singlets. However, if = (1,2), then a = b and 
the part in Eq. (18) which is proportional to (a — b)/2 vanishes. 

If in the combination |0i| 2 |0j| 2 in Eq. (13) the fields (<f>i,(j>j) form one of the 
pairs (0i, 02), (03, 04), or (05,06), then the combination |0«| 2 |0j| 2 is a total 
singlet (on its own) and it can be written directly into the scalar potential as 
a|0j| 2 |0j| 2 , where a is some real- valued constant. 

Moreover, the symmetries V 2 and V 3 enforce the products |0j| 2 and |0j| 4 [in 
Eq. (13)] to appear in the scalar potential only as 



A (l0il 2 + N 2 ) + A (l0sl 2 + N 2 ) + & (I0 5 | 2 + \M 2 ) 

+ «1 (|0l| 4 + |0 2 | 4 ) + «2 (|03| 4 + |04| 4 ) + *3 (|0 5 | 4 + |0 6 | 4 ) , (19) 

where /if, [/%, Ki, k 2 , and k 3 are real-valued constants. In total, the most 
general scalar potential involving only the unprimed fields, but neither the 
fields 07, 08, 09, and 0io nor the product \9\ 2 , reads 



Vx=fi (I0il 2 + 10 2 | 2 ) + 14 (I0 3 | 2 + 10 4 | 2 ) + 14 (I0 5 | 2 + 10 6 | 2 ) 



+ «1 (|01 



+ 



+ 



+ 



+ rf 3 
+ 4 



+ |04| 2 ) (|0 5 | 2 +|0 6 | 2 ) + 



1^2 | 



3^4 



4 1 

d 2 

2 + d; 



+ «3 
0l| 2 - 



+ 



106^ 



- I06| 2 ) 



+ d 7 |0 5 06 

+d 9 [(0|0 4 ) 2 + (0i0 3 ) 2 ] + d 10 [(0^0 6 ) 2 + (0^05) 

+ d n (0l0 2 + 2 0l) (0^04 + 0403 ) 
+ d 12 (0l0 2 + 0^0x) (050 6 + 6 5 ) 
+ di 3 (0^04 + 040 3 ) (0506 + 06^) 
+ d U (0304 - 0403) (0506 ~ 060s) , 



(0l0 2 ) 2 + (020l) S 



(20) 



where di,d 2 , . . . ,du are real- valued constants. We will assume that ki,k 2 ,k 3 > 
and dn,di3 > and we will choose all other coupling constants to be neg- 
ative. Again, we observe that the first nine terms in Eq. (20) exhibit three 
accidental U(l) symmetries, which act on the pairs of VEVs ((0i), (02)), 
((03), (04)), an d ((05), (06)), respectively. The rest of the terms in the po- 
tential V\ break these symmetries and will therefore determine the vacuum 
alignment mechanism of the fields. First, we note that the term with the co- 
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efficient <i 4 tends (for large values of \d\\) to induce a splitting between \((j>z)\ 
and |(04)| as well as between |(0s)| and |(06)|- Second, we observe that the 
term with the coefficient d u has (for large values of |d 14 |) the tendency to 
trigger relative phases (different from and ir) between (0 3 ) and (0 4 ) as well 
as between (0 5 ) and (06)- However, if we require that 

d 6 d 7 > Ad\ and d^di® > 4g^ 4 , (21) 

then the potential V\ is minimized by the VEVs of the unprimed fields, which 
are pairwise degenerate in their magnitudes, i.e., they satisfy 

l(0i)l = l(0 2 )l, 1(03)1 = 1(04)1, |(0 5 )| = |(0 6 )| (22a) 
and are also pairwise relatively real, i.e., 

(0l) (03) (05) 



(02)' (04)' (06 



e{-l,l}, (22b) 



where, in addition, the choice d\i < and dn,di 3 > implies a correlation 
between the different pairs of VEVs in terms of 

(0 2 )(06) _1 ^ (0 2 )(04)" (22C) 

i.e., the relative sign between (0i) and (02) is equal to the relative sign between 
(05) and (06) and opposite to the relative sign between (0 3 ) and (0 4 ). 



3. 6 The potential of the fields 0' 1; (f>' 2 , . . . , 0' 6 



The only two-fold products of the primed fields 4>[ (« = 1, 2, . . . , 6) which are 
allowed by the discrete symmetries Vi,V 5 , and V e are the absolute squares 
|0£| 2 . Furthermore, the permutation symmetries V 2 and V 3 yield for the most 
general dimension-two terms of the primed fields the expression 

m\ 2 + i0 2 i 2 ) + m\ 2 + i0 4 i 2 ) + i4m\ 2 + w,\% (23) 

where n\, n\, and /Zg are real-valued constants. The permutation symmetries 
V 2 and V 3 additionally require the most general products of the absolute 
squares |0-| 2 to be 



^4 (10'il 2 + |0' 2 | 2 ) 2 + « 5 (|0 3 | 2 + |0 4 | 2 ) 2 + ^ (|0 5 | 2 + |0 6 | 2 ) 2 

+«i (m 2 + i0' 2 i 2 ) (m 2 + m 2 ) + 02 (i0;i 2 + 10'/) (105I 2 + w) 
+« 3 (i0 3 r + i0 4 i 2 ) (i0 5 i 2 + w) + 04 (103I 2 - i0 4 i 2 ) (i0 5 i 2 - m 2 ) 

+ a 5 |#ty 2 | 2 + a6|0 3 V 4 | 2 + a 7 \<f>' 5 UX , (24) 
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where K4, k§, kq, and a±, a 2 , . . . ,07 are real-valued constants. Note that three- 
fold products of the primed fields are forbidden by the £V c harge assignment. 
The discrete symmetry D 6 requires that the remaining interactions of the 
primed fields can be written as products of the operators 4>'\4> 2 , 4>'\4>'a^ an d 
0'50'g (and their complex conjugates). Hence, the operators which involve the 
fields of only one of the pairs ((/>[, <f>' 2 ), (0' 3 , 4 ), or (0' 5 , 0' 6 ) are restricted by the 
symmetries T> 3 and V4 to be on the form 



a 8 
+ ai 



(0;V 2 ) 2 + (0 2 V'i) : 



/t. 



a 9 



2a 8 K (<P'W 2 ) 2 + 2a 9 K (0' 3 V 4 ) 2 + 2a 10 K (0' 5 T 0' 6 ) 



4) 2 



(#V 3 ) S 



n't. 



(25) 



where a 8 , ag, and aio are real- valued constants. Furthermore, the symmetries 
V2,T> 3 ,V 5 , and V 6 restrict the only dimension-four terms involving at least 
three different primed fields (or their complex conjugates) to be 



a n (0; V 2 + 2 f 0i) (0 3 f 4 + 04 Vs) 
+ a 12 (0; V' 2 + 2 Vi) (05 Ve + ^6 V 6 ) 
+ ai3 (^sVl + 01*03 ) (05*06 + $$V 5 ) 

+ ai4 (03 V4 - 04 V3) (05 f 06 - 06 t 05) 

= 4a 11 K(0' 1 t 2 M03 t 04) + 4a 12 3?(0 / 1 t 2 )^(0 5 t , 6 ) 

+ 4a 13 K(03 t 04)^(05 t , 6) - 4a 14 3(0 3 V4)3(05 t , 6), (26) 



where an, ai 2 , 013, and 014 are real- valued constants. Taking everything into 
account, the most general scalar potential involving only the primed fields <$>\ 
reads 
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'el , 



vi=ii 2 M\ 2 + \m + dm 2 + 1011 2 ) + m\ 2 + 10 

+ k 4 (|0;i 2 + 10' 2 | 2 ) 2 + K5 (|0' 3 | 2 + |0;i 2 ) 2 + K6 (|0' 5 | 2 + 
+01(10; 

+ a 5 |0' 1 t 0' 2 | 2 + a 6 |^V4| 2 + arl^'sVel 2 + a 8 [(0'iV 2 ) 2 + (<t>'U'if 
+ a 9 [(03 V 4 ) 2 + (04V3) 2 ] + aio [(05 t 6 ) 2 + (0 , 6 t 0' 5 ) 5 
+ an 
+ a i2 



0;i 2 + i0 2 i 2 ) 2 + «5(i^i 2 + i^i 2 ) 2 + **m\ 2 + i0ei 2 ) 2 

^il 2 + I0 2 | 2 )(I0' 3 | 2 + \m + a 2 (|0;| 2 + |0' 2 | 2 ) (|0' 5 | 2 + |0' 6 | 2 ) 
10'sl 2 + I0 4 | 2 ) (I0' 5 | 2 + W) + «4 (|0' 3 | 2 - I0' 4 | 2 ) (I0' 5 | 2 - I0 6 l 
'M\ 2 + a^U'f + *rWM\ 2 + [(<UM) 2 + (^V'i) 2 



y/02 + 2 Vl) (03 t 0l + 4 V 3/ 
(0^02 + 2 t 0'l) (0^0; + 6 V 5 ) 



+ A13 (0^04 + V 3 ) (0^06 + 6 V E 



+ ai4 



(0^0; - 0IV3) (0^06 - 0^05) • (27) 

In Eq. (27), we will assume that k&,k$,kq > and we will choose all other 
coupling constants to be negative. As in the discussion of the potential Vi, 
we observe that the first nine terms in Eq. (27) exhibit three accidental U(l) 
symmetries, which act on the pairs of VEVs ((0'i), (0' 2 )), ((03), (04)), an d 
((0' 5 ), (0g)), respectively. The rest of the terms in the potential V{ break these 
symmetries and will therefore determine the vacuum alignment mechanism of 
the fields. If we, in analogy to the potential Vi, require that 

a 6 aj > 4a 4 and a 9 a w > Aa 1A , (28) 

then the potential V[ is minimized by the VEVs of the primed fields, which 
are pairwise degenerate in their magnitudes, i.e., they satisfy 

l(0'i)l = l(0 2 )l, 1(03)1 = 1(01)1, 1(05)1 = 1(06)1 (29a) 
and are also pairwise relatively real obeying 

(0'i) _ (0' 3 ) _ (0' 5 ) 



(0' 2 ) (00 (0' 6 ) 



e{-l,l}. (29b) 



Note that in Eq. (29b) the pairs of the VEVs ((0;), (0' 2 )), ((0 3 ), (0^)), and 
((0' 5 ), (0g)) have the same relative phase, i.e., the VEVs in the pairs are either 
all oriented parallel or all oriented antiparallel. 



3. 7 Mixing among the fields 0' 1; 0' 2 , ... , 6 and 



The discrete symmetry D 6 requires all renormalizable terms mixing the primed 
fields 0^ or 0^ (i = 1, 2, . . . , 6) with the unprimed fields 0-,- or 0] (j = 
1, 2, . . . , 6) to have an even mass dimension. Taking the combinations in 
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Eq. (7) and the product 4>\ 4 03,4 into account (which are all £> 6 -singlets), we 
obtain that the Dx-invariant operator products, which mix the fields $\ and 
(f>j, must be of the types 

01^01,201,2, K 0^3,403,4 , 01^05,605,6, (30) 

where i±, i 2 — 1, 2, . . . , 6. The symmetry P 4 , which acts only on the unprimed 
fields, requires the combinations in Eq. (30) to be on the form 0l 1 0^ 2 |0j| 2 , 
where j = 1,2,... ,6. Next, the symmetries V 5 and V 6 imply that the opera- 
tors in Eq. (30) are all in fact ~ |0'i| 2 |0:/| 2 - As a result, the most general renor- 
malizable interactions of the fields </>[, (f>' 2 , . . . , 0g with the fields 0i, 02, . . . ,06 
are 



^ 2 =(|0 / 1 | 2 + |0 2 | 2 ) 

x h (|0 X | 2 + |0 2 | 2 ) + b 2 (|0 3 | 2 + |0 4 | 2 ) + h (|0 5 | 2 + |0 6 | 2 ) 

+ & 4 (|0 , 1 | 2 -|0 2 | 2 ) (|0i| 2 -|0 2 | 2 ) 

+ (|0 3 | 2 + I0 4 l 



& 5 (|0l| 2 + 



+ (|0' 3 | 2 " I0 4 | 2 ) [68 (10 
+ (|05I 2 + I0 6 | 2 ) 
1 



+ i0 4 r)+& 7 (i05i 2 +i0 6 i 2 ) 
- 4 r)+6 9 (i0 5 i 2 -i0 6 i 2 ) 



X 



'lO 



+ |0 2 | 2 ) + 6ll (|0 3 | 2 + |04| 2 ) + 6l2 (|0 5 | 2 + 



+ (I05I 2 " 



/ |2 



'13 



|0 4 ] 2 )+6l4(|0 5 | 2 -|06 



(31) 



where bi,b 2 , . . . , &14 are real- valued constants. In Eq. (31), we will assume all 
coupling constants to be negative. 

In order to recover the (same) vacuum alignment mechanism that is operative 
for the potentials V 1: V[, and V 7j8 also for the full SM singlet scalar potential 
V = V\ + V[ + V 2 + Vt-,8, we will have to ensure that the mixed terms in 
the potential V 2 do not induce a splitting between the pairwise degenerate 
magnitudes of the VEVs. If we require the coupling constants in the potentials 
Vy, V{, and V 2 to fulfill 

a 5 d 5 > 464, clqcIq > 46g, a 6 d 7 > 46g, a 7 d 6 > Abf 3 , a 7 d 7 > Ab\ 4 , (32) 

then the total multi-scalar potential V is indeed minimized by the VEVs of 
Eqs. (12), (22), and (29). 

We will suppose that all of the SM singlet scalar fields break the flavor sym- 
metries by acquiring their VEVs at a high mass scale (somewhat below the 
fundamental mass scale Mi), and thereby, giving rise to a small expansion 
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parameter 



where % = 1,2,... ,10 and j = 1,2,... ,6. Such small hierarchies can arise 
from large hierarchies in supersymmetric theories when the scalar fields acquire 
their VEVs along a "D-flat" direction [67,68]. 

In the next two sections, Sees. 4 and 5, we will discuss the (effective) Yukawa 
interactions of the leptons that will lead to the mass matrix textures of the 
leptons, and eventually, to the masses and mixings of them, which will be 
discussed in Sees. 6 and 7. 



4 Yukawa interactions of the charged leptons 



Consider the effective Yukawa coupling operators 0^ which generate the en- 
tries in the charged lepton mass matrix via the mass terms 

S?£ = L;H 2 e afi E l3 + h.c., (34) 

where a, (3 = e, \i, r. (See Fig. 2.) We will denote the total number of times that 




Fig. 2. Non-renormalizable terms generating the effective Yukawa couplings in the 
matrix (O e a/3 ). 

the fields <f>i, 2 , ■ ■ ■ ,06 appear in the operator by n\ and the total number 
of times that their complex conjugates <f>\, (f) 2 , • • ■ ,06 appear in the operator 
@uf3 by n 2 . Now, invariance under transformation of the discrete symmetry 
T>i implies that for the first column of the Yukawa interaction matrix (^U), 
i.e., for /3 — e, it must hold that ri\ — n 2 = 4. For the second and third column 
of the Yukawa interaction matrix, i.e., for (5 = fj,,T, the discrete symmetry 
T>\ instead requires that n\ — n 2 = 1. In addition, we conclude from the 
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transformation properties of the fundamental Froggatt-Nielsen states under 
the discrete symmetry D 6 that the operators ^ and <^ T , where a = e,fi,T, 
can neither involve the field 0g nor the field 0io- This is, however, not true 
for the operators G l ae (a = e, /i, r) in the first column of the effective Yukawa 
coupling matrix. 



4-1 The first row and column of the charged lepton mass matrix 



Invariance under transformations of the U(l) symmetries requires the U(l) 
charges of the entries 6* an<4 - m the ^ TS ^ row °f the effective Yukawa 
coupling matrix {& a a) to be (1, —1,0). Since the fields 0g and 0io cannot be 
involved in the generation of the e-fi- and e-r-elements of the charged lepton 
mass matrix, the U(l) charge assignment immediately implies that any mass 
operator giving rise to these e-/i- and e-r-elements must involve the term 
~ 0i,2# 2 /(^i) 3 - Next, the symmetries V 5 and V 6 yield to leading order for 
the operators & and & er the two possible terms ~ 0i0' 1 # 2 /(Mi) 4 and ~ 
0202^ 2 /(^i) 4 - In conjuction with the requirement ri\ — ri2 = 1, the symmetry 
Vq implies that any further operators contributing to the operator tf* or 
G l eT must have at least two powers of mass dimension more than the terms 
0i0^ 2 /( M i) 4 and 0202#7( M i) 4 - We wil1 therefore neglect these additional 
operators. 

From the transformation properties of the right-handed electron E e and the 
fundamental Froggatt-Nielsen states under transformations of the discrete 
symmetries V\ and T> 6 , we conclude that the operators 6 l ee , , and @ l re in 
the first column of the effective Yukawa coupling matrix must involve at least 
a four-fold product of fields taken from the set {0i, 02, 03, 04, 05, 4>e} times a 
field taken from the set {0 9 ,0i O }. Possible lowest-dimensional contributions 
to the operator ff e ee , which are consistent with the symmetries of our model, 
are, e.g., given by ~ 1O [(0 3 ) 4 + (0 4 ) 4 ] /(Mi) 5 and ~ 1O (0 3 ) 2 (0 4 )7(M 1 ) 5 . For 
brevity, we will take the operator 

^ 010 



(0s) 4 + (0 4 ) 4 (35) 



as a representative of these contributions. The remaining operators & and 



[Mi 



& Te have a mass dimension that is greater than or equal to the mass dimension 
of the terms in Eq. (35). However, the effects of these terms on the leptonic 
mixing angles will turn out to be negligible in comparison with the contribu- 
tions coming from other entries of the charged lepton mass matrix. 

In total, the first row of the effective Yukawa coupling matrix of the charged 
leptons, which is consistent with all of the discrete symmetries, is to leading 
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order 

«J = (a y [(0 3 ) 4 + (0 4 ) 4 ] B 1 - 020' 2 ] B 1 [M'l + <h<&]) ■ ( 36a ) 
Here the dimensionful coefficients A 1 and B\ are given by 



* =1 W < 36b > 
B >= y W < 36c > 

where the quantities Y„ and are arbitrary order unity coefficients and Mi 
is the high mass scale of the intermediate Froggatt-Nielsen states. Note that 
at the level of the fundamental theory, the permutation symmetry V4, which 
interchanges the second and third generations of the leptons, is also propagated 
to the heavy Froggatt-Nielsen states. This establishes a degeneracy of the 
associated Yukawa couplings and the explicit masses of these states, which 
is then translated into a degeneracy of the corresponding effective Yukawa 
couplings of the low-energy theory. 



4-2 The 2-3-submatrix of the charged lepton mass matrix 



In the 2-3-submatrix of the charged lepton mass matrix, the U(l) charges 
of the operators { a ->fl = A 4 5 r ) m ust be (0,0,0). The lowest dimensional 
operators which fulfill this condition as well as the constraint n\ — n 2 = 1 
are proportional to 03,4/Mi or 05 j6 #/(Mi) 2 . Furthermore, invariance under 
transformation of the discrete symmetries D 5 and T>$ implies that the lowest 
dimensional operators 6^ in the 2-3-submatrix with n\ — n 2 = 1 are of the 
types 



0404 05050 0606^ ^ 



(Mx) 2 ' (Mi) 2 ' (Mx)3' (Mx)3 

Thus, the most general 2-3-submatrix of the matrix which involves only 

these combinations and is invariant under transformations of the remaining 
discrete symmetries, is found to be 

/ C(0' 3 3 - 0' 4 4 ) + £(0505 - 0' 6 6 ) 

C<(0303 + 4 4 ) +^(0505 + 6 6 ), 

(38a) 

Here the dimensionful coefficients C and D are given by 
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c=5 w (38b) 
D=y w < 38c) 

where the quantities Y s c and are arbitrary order unity coefficients. Actually, 
contributions to the next-leading operators 6* and O^. are, e.g., given by 



H-t: 
T-fi: 



( M iY 



(^10708 + O20807) (03 + 04)03, (39a) 



(Mi, 



^(^20708 + ai0 8 7 )(03 ~ <Pi)<P'^ ( 39b ) 



where a\ and 02 are some complex-valued constants. Note that the terms in 
Eqs. (39) carry only one unit of mass dimension more than, e.g., the term 
/ 5 5 6'/(M 1 ) 3 in Eq. (38a). When diagonalizing the mass matrix, however, it 
will turn out that the associated corrections to the leptonic mixing parameters 
and the lepton masses are in fact negligible. 



4-3 The charged lepton mass matrix 



Combining the results of Sees. 4.1 and 4.2, the leading order effective Yukawa 
coupling matrix of the charged leptons is 



[(4> 3 f + W4) 4 ] B 1 - <fe<#, 2 ] Si + 4>2<P' 2 ] 

C(03*3 - ¥i<t>i) + D(<l>' 5 <t>5 - <t>' e 4>6) 

C(4>' 3 4> 3 + 0^4) + B(<#>5<fe 



(40) 



where the dimensionful couplings A 1 , B x , C, and D are given in Eqs. (36b), 
(36c), (38b), and (38c), respectively. Inserting the VEVs in Eqs. (29) and (22) 
into the corresponding operators of the matrix (^£3), we observe that, due to 
the vacuum alignment mechanism of the SM singlet scalar fields, in some of 
the entries of the matrix {{M() a p), the spontaneously generated effective mass 
terms of a given order exactly cancel, whereas in other entries they do not. 
Furthermore, the vacuum alignment mechanism correlates these cancellations 
in the different entries of the matrix ((M^) a/ g) in such a way that after SSB the 
charged lepton mass matrix Mg can be of the two possible asymmetric forms 



m 



cxp 



e 3 e 2 e 4 ^ 



e 3 e e 2 



e 3 e 2 1 



(41) 
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and 



/ e 3 £ 4 £ 2\ 



e 3 1 e 2 



e 3 e 2 e 



(42) 



where we have also introduced the appropriate orders of magnitude for the 
matrix elements (M^)^ e and (Mg) Te as well as for the "phenomenological" (in 
contrast to exact) texture zeros arising in Eqs. (36a) and (38a). Here m^. xp is 
the (experimental) tau mass. Note that a permutation of the second and third 
generations L M <-> L T , <-> _E r leads from one solution to another. Let us 
consider the first one for our remaining discussion. 



5 Yukawa interactions of the neutrinos 



Consider the effective Yukawa coupling operators G v a ^ which generate the en- 
tries in the neutrino mass matrix M v via the mass terms 

^y=^C^ + h.c., (43) 

where a, (3 = e, fi,r and M 2 is the relevant high mass scale which is responsible 
for the smallness of the neutrino masses in comparison with the charged lepton 
masses. Since the SM singlet neutrinos as well as the Higgs doublets are T>\- 
singlets, the presence of the fields 0i, 02, • • • ,06 (which transform non-trivially 
under the discrete symmetry Vf) in the operators 0^$ * s forbidden. Hence, 
the only scalar fields that can be involved in the leading order operators 
are 7 , 8 , 9 , 0i O , and 9. 



5. 1 Effective Yukawa interactions of the neutrinos 



The operators 0^ T , , and 6" T must have the U(l) charge structure 
(0, —2, 0). An example of the lowest dimensional operators which achieve this 



is 



0901Q0 

(M x ) 5 



(4) 2 + (4f 



(44) 
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Since the U(l) charge structure of the entry ff v ee is (—2, 0, 0), the operator G v ee 
is to leading order 



Y 



a (M ± y 







(45) 



where Y£ is an arbitrary order unity coefficient. Note that the operators of the 
type in Eq. (44) carry two units of mass dimension more than the operator 



ee' 



The operators and G v er (as well as G v and <^ e ) must have the U(l) charge 
structure (—1,-1,0). Therefore, the lowest dimensional terms which con- 
tribute to these operators are ~ fa /Mi and ~ fa /Mi. The symmetries V 2 , ^3, 
and V 4 then yield to leading order ^ M = Y b u fa/Mi and G v er = Y h v fa/M u 
where Y b u is an arbitrary order unity coefficient. Note that the operator G v ee 
in Eq. (45) carries two units of mass dimension more than the operators 
&en ~ 4>7/Mi and G v eT ~ fa/ Mi. We will therefore not in detail consider 
the structure of the highly suppressed terms that appear in the 2-3-submatrix 
of the neutrino mass matrix as the operator given in Eq. (44). 

In total, the most general effective Yukawa coupling matrix of the neutrinos 
i^ap) tnat * s consistent with the symmetries of our model is to leading order 
given by 



^ A 2 B 2 



O B 2 
B 3 

Here the dimensionful coefficients A 2 , B 2 , and B 3 are given by 



(46) 



/ 



fafao9 

A2 ~ Ya (MO 3 ' { } 

8* = *?^, (47b) 

Bs = Y b ^, (47c) 

where the quantity Y£ is an order unity coefficient. The leading order tree level 
realizations of the higher- dimensional operators, which generate the neutrino 
masses, are shown in Figs. 3 and 4. Note that the coefficients B 2 and _B 3 
contain the same Yukawa coupling constant Y^. Furthermore, we point out 
that the texture zeros in the 2-3-submatrix of the effective neutrino Yukawa 
matrix should be understood as "phenomenological" ones, since they actually 
represent highly suppressed operators carrying two units of mass dimension 
more than the entry ff v ee . 



23 



1 

(f>a 



Fig. 3. The dimension six operator for a = (j,,t and 4>^ = <pj, 4> T = 4>8 generating 
the e-(i- and e-r-elements in the effective neutrino mass matrix. 



Ff 



II, 



4>io 



iff! 



Ff 



II, 



II, 



Fig. 4. The dimension eight operators generating the e-e-element in the effective 
neutrino mass matrix. 

5.2 The neutrino mass matrix 



Inserting the VEVs in Eq. (12) into the effective operators in Eqs. (47), we 
obtain from Eq. (46) the effective neutrino mass matrix (to leading order) 



Mo 



1 Y:e> Y^e 



Yfe e 
±Y h u e e 



(48) 



where we have introduced the actual orders of magnitude of the higher-order 
corrections to the texture zeros in the 2-3-submatrix of the matrix in Eq. (46). 
Note that after SSB the symmetries determine the e-pL- and e-r-elements to be 
exactly degenerate (up to a sign), giving rise to an atmospheric mixing angle 
which is close to maximal (higher-order corrections to exact maximal atmo- 
spheric mixing mainly come from the fi-r- and r-/i-elements of the charged 
lepton mass matrix). Introducing an "absolute" neutrino mass scale m u and 



24 



choosing Y£ /Y£ ~ 1, we can write the neutrino mass matrix in Eq. (48) as 



(49) 



where 

™. = ^K, (50) 

Note that we have chosen the minus signs for the e-r- and r-e-elements due 
to our freedom of absorbing the corresponding phase into the order unity 
coefficients in the charged lepton sector. Furthermore, it is important to keep 
in mind that the entries "1" and "—1" of the matrix in Eq. (49) indeed denote 
matrix elements, which are degenerate to a very high accuracy, whereas the 
other entries are only known up to their order unity coefficients. 



s z 1 -1 



1 e 4 e 4 
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6 Lepton masses and leptonic mixings 



From the results of the last two sections we have seen that the lepton mass 
matrices are given by 



M, 



i - m T 



cxp 



/ ,3 £ 2 e 4\ 



e 3 e e 2 



e A e 2 1 



/ 



and M v ~ m v 



1 -1 



1 4 

lee 
e 



1 e 4 ^ 



where again m® xp is the (experimental) tau mass, m u is some "absolute" neu- 
trino mass scale, e ~ 0.1 is the small expansion parameter (defined in Sec. 3), 
and only the order of magnitude of the matrix elements have been indicated. 
In the above expression, the first matrix is the charged lepton mass matrix 
and the second matrix is the neutrino mass matrix. Note that the small expan- 
sion parameter e is the same for both matrices. In order to find the leptonic 
mixing angles and the lepton masses, we have to perform diagonalizations of 
the above two displayed mass matrices. Let us begin with the diagonalization 
of the charged lepton mass matrix. Since this matrix is not symmetric and we 
want to extract the relevant mass and mixing parameters, we have, in fact, 
to diagonalize the matrix product MiM\, which is a symmetric matrix. Bi- 
unitary diagonalization of the matrix Me implies that \j\MgVi = Ale, where 
Aie = diag (m e ,m M ,m T ) is the diagonalized charged lepton mass matrix con- 
taining the masses of the charged leptons (i.e., the eigenvalues of the matrix 
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\J M s M\~) and Ui and Vt are two unitary mixing matrices. Thus, we obtain 

M.iM\ = U\~MiiM\1) r £, where Ug will be the charged lepton mixing matrix. 
Note that the mixing matrix Ve does not appear in the diagonalization of the 
matrix product M(M[. Actually, since the matrix Mi is a real matrix, we have 
MiM\ = MeMj = Uj M t MjUi. Next, straight-forward calculations (using 
App. A) yield 

M e M\ = (m^) 2 diag (e 6 - 2e 7 + 0(e 9 ), e 2 + e 4 + G(e 5 ), 1 + 2e 4 + 0(e 5 )) 

(51) 



and 



U e = {{U t ) aa ) ~ 



-0.995 -0.0997 -0.000212 
0.0997 -0.995 -0.0111 
-0.000897 0.0111 -1.000 



(52) 



such that UjUe ~ UeUj ~ I3, where a — e, /x, r and a = 1, 2, 3. Note that the 
charged lepton mixing matrix is independent of the absolute charged lepton 
mass scale m^. xp , i.e., the tau mass. This means that the charged lepton (and, 
of course, the leptonic) mixing angles will not depend on the tau mass. The 
charged lepton masses are given as the square roots of (the absolute values of) 
the eigenvalues of the matrix product M(M t 



m, 



rn 



cxp 



is of the matrix product M^Mj, i.e., m e = m^e 3 (1 — e + 0(e 2 )), 
(l + \e 2 + C(e 3 )J, and m T = rrf^ (1 + e 4 + C(e 5 )). Thus, the 
order-of-magnitude relations for the charged lepton masses are given by 



m e /m M ~ e 2 ~ 10 2 and m^/m T ~ e ~ 10 



(53) 



li J cxp — 



which approximately fit the experimentally observed values, i. e., (m e /m 
4.8 ■ 10 -3 and (m M /m r ) cxp ~ 5.9 ■ 10~ 2 [28]. Furthermore, the charged lepton 
mixing angles (in the "standard" parameterization) are found to be (using 
App. A) [69,70,49] 



6 12 = arctan 



i)e2 



(U e ) el 
6>f 3 = arcsm|(L^) e3 | 



6 23 = arctan 



t3 



--6-^ + 0(6% 

2e 2 + e 5 + C(e 6 ), 
= e 2 + e 3 + 0( e 4 ). 



(54) 
(55) 
(56) 



'12 



5.72 c 



6°, 9[ 
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Inserting e ~ 0.1 into Eqs. (54) - (56), we obtain 
0.0122° 0.01°, and ^3 — 0.637° ~ 0.6°, i.e., the mixing angles in the charged 
lepton sector are all small. The mixing angle 9\ 2 is the only one that is not 
negligible. Thus, it will finally yield a non-zero contribution to the leptonic 
mixing angles [49]. 
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Next, we want to diagonalize the neutrino mass matrix M v . The diagonal- 
ization procedure of the matrix M v is easier than the one of the matrix Mi, 
since the matrix M v is a symmetric matrix. Diagonalization of the matrix M v 
directly implies that UjM u U u = M. v , where M. u = diag(m Ul ,m U2 ,m U3 ) is 
the diagonalized neutrino mass matrix containing the masses of the neutrino 
mass eigenstates (i.e., the eigenvalues of the matrix M v ) and U v is the unitary 
neutrino mixing matrix. Straight-forward calculations (using App. A) yield 



M„ = m u diag Q (e 2 + V8T^) }- (e 2 - + , 2e l 
= m v diag (^/2 + ^e 2 + ^=e 4 + 0(e 8 ),-V2 + y 



8V% 



e 4 + £>(e 8 ),2e 4 



and 



U v = ((U v ) aa ) 



( 0.708 0.706 ^ 

0.499 -0.501 ^= 
v -0.499 0.501 ^, 



(58) 



such that U V U V ~ U v Ul ~ 1 3 , where a = e,/j,,r and a = 1,2,3. Note that 
the neutrino mixing matrix U u is independent of the absolute neutrino mass 
scale m v . This means that the neutrino (and, of course, the leptonic) mixing 
angles will not depend on this scale. The (physical) neutrino masses are given 
as (the absolute values of) the eigenvalues of the neutrino mass matrix M v , 
i.e., mi = \m Vl \ ~ \[2m v , m 2 = \m U2 \ ~ \plm v , and m 3 = \m V3 \ = 2e 4 m l/ ss 0, 
which means that we have an inverted neutrino mass hierarchy spectrum (i.e., 
m 3 <C m x ~ m 2 =>- ~ (Am 2 ,! <C |Am§ 2 | ^ jAm 2 ,!). Naturally, the 



neutrino mass squared differences Am 2 5 = m 2 a — m 2 ft = m 2 — m\, where m Va 
is the mass and m a is the physical mass of the ath neutrino mass eigenstate, 
respectively, are given as follows: 



(57) 



Am 2 , = ml (-2V2e 2 - + 0(e 10 )^ , (59) 

Am 2 2 = ml (-2 + V2e 2 - ^ + 0(e 6 )) , (60) 
Am 2 , = ml (-2 - v^e 2 - ^ + G(e 6 )) . (61) 

Note that the leading order terms in Eqs. (59) - (61) are all negative, which is 
natural, since we have an inverted mass hierarchy spectrum for the neutrinos. 
Hence, the solar and atmospheric neutrino mass squared differences are given 
by Am| = (Am 2 ,! ~ 2y / 2e 2 m 2 , and Am 2 tm = |Am| 2 | ~ (Am 2 ,) ~ 2m 2 ,, 
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respectively. The "experimental" values of these quantities are Am^ ~ 5.0 • 
KT 5 eV 2 [23] and Am 2 tm ~ 2.5 • 10" 3 eV 2 [81-PT1 Thus, we extract that m v ~ 
0.04 eV and e ~ 0.1, which is consistent and agrees perfectly with our used 
value for the small expansion parameter e. In other words, using e ~ 0.1 and 
m u ~ 0.04 eV for the absolute neutrino mass scale (Nobody knows where this 
value comes from!), we obtain the presently preferred values for the solar and 
atmospheric neutrino mass squared differences. Furthermore, we have that 
rrii — m 2 — 0.05 eV and m 3 ~ 1 ■ 10~ 5 eV ~ 0. Similarly, as for the charged 
leptons, the neutrino mixing angles (in the "standard" parameterization) are 
found to be (using App. A) [69,70,49] 



d\ 2 = arctan 



(U, 



v e2 



v el 



6>i3 = arcsin \{U v )es\ 
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: arctan 



{U V ) TZ 



7T 1 

4 + 472 f 



7T 

4" 



96^ 



+ 0{e 



10\ 



(62) 
(63) 
(64) 



Inserting e ~ 0.1 into Eqs. (62) - (64), we obtain 6 v l2 ~ 45.1° w 45°, 9% 3 = 0, 
and ^23 = 45°, which means that we have (nearly) bimaximal mixing in the 
neutrino sector. Note that the values for Q\ z and #23 are exact, since they are 
only determined from the matrix elements of the third column of the matrix 
U v in Eq. (58). 



7 The leptonic mixing angles 



The leptonic mixing angles (or parameters if one also includes the CP violation 
phase 5) are given by the leptonic mixing matrix^]. The leptonic mixing matrix 
U is composed of the charged lepton mixing matrix Ue and the neutrino mixing 
matrix U v as follows: 

U = U\U U . (65) 

The matrix Ui (U v ) rotates the left-handed charged lepton fields (the neutrino 
fields) so that the charged lepton mass matrix Mi (the neutrino mass matrix 
My) becomes diagonal (see Sec. 6), i.e., it relates the flavor state and mass 
eigenstate bases. Thus, one can look upon the matrix U [U a b = (UjU u ) a b = 

6 @ 99.73% C.L.: 2.3- 1(T 5 eV 2 < Am| < 3.7- l(T 4 eV 2 [23]; best-fit: Am| ~ 
5.0 • 10~ 5 eV 2 [23] 

7 @ 90% C.L.: 1.6 • 10" 3 eV 2 < Am^ < 3.9 • 10" 3 eV 2 [9]; best-fit: Am 2 tni ~ 
2.5 • 10" 3 eV 2 [8] 

° The leptonic mixing matrix is sometimes called the Maki-Nakag awa S akat a 
(MNS) mixing matrix [71]. 
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Ea=e,nA U })aa(U u ) ab = E a =e, M ,r(^)aa(^)a&] as the messenger between the 
mass eigenstate bases of the charged leptons and the neutrinos. Inserting the 
matrices in Eqs. (52) and (58) into the definition in Eq. (65), we find that 



^ -0.655 -0.753 0.0699^ 



U = (U ab ) 



-0.573 0.434 -0.696 
0.493 -0.495 -0.715 



(66) 



which fulfill the condition WU = UW — I3 to a very good accuracy. Now, 
the leptonic mixing angles (in the "standard" parameterization [see App. A]) 
can be read off as follows [69,70,49]: 



U 



12 



912 = arctan 

9 13 = arcsin \ Ui 3 \, 
9 2 s = arctan 
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(67) 
(68) 
(69) 



Thus, inserting the appropriate matrix elements of the matrix U expressed in 
terms of the small expansion parameter e (e <C 1) in Eq. (65) into Eqs. (67) - 
(69), we obtain 



e «=T2-^r/ + ° {£l) ' (71) 

fe = f-^-f' + 0( e '). (72) 

Note that the first correction to the atmospheric (neutrino) mixing angle #23 
is of second order in the small expansion parameter e, and it is therefore 
very small, i.e., the atmospheric mixing angle stays nearly maximal, #23 — f ■ 
However, the first corrections to the solar (neutrino) mixing angle #12 and the 
reactor (neutrino) mixing angle (the so-called CHOOZ mixing angle) #13 are 
both of first order in the small expansion parameter e and they are of exactly 
the same size, but with opposite sign. Thus, we have the first order relation 
9 12 c± I — #i 3 . Finally, inserting e ~ 0.1 in Eqs. (70) - (72), we find that the 
leptonic mixing angles are 

012 41.0° w 41°, 13 ~ 4.01° w 4°, and 6 23 ~ 44.2° w 44°, 

which means that we have bilarge leptonic mixing. These values of the leptonic 
mixing angles lie within the ranges preferred by the MSW LMA solution of the 
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solar neutrino problem^, atmospheric neutrino dataf 10 ] (nearly maximal atmo- 
spheric mixing), and CHOOZ reactor neutrino data. The so-called CHOOZ 
upper bound is sin 2 6 Vi < 0.10 {i.e., \6 Vi \ < 9.2°) [72-74]. In particular, the 
obtained value for the solar mixing angle # 12 implies a significant deviation 
from maximal solar mixing. However, the solar mixing angle is bounded from 
below by approximately 41° and it is therefore still too close to maximal to 
be in the 90% (or 95% or 99%) confidence level region of the MSW LMA 
solution [23]. 

As we have seen the symmetries of our model generate the hierarchical charged 
lepton mass spectrum as well as an essentially maximal atmospheric mixing 
angle. However, these symmetries only determine the order of magnitude of 
the entries (Mf) aj 3 in the charged lepton mass matrix M#. In order to decide 
whether the model can naturally give the MSW LMA solution or not, we 
can test the robustness of the above calculated leptonic mixing angles (and 
charged lepton masses) under variation of the involved order unity coefficients. 
For example, by changing the ratio of the order unity coefficients Y^/Yj from 
1 to 2 leads to 

9 12 ~37°, 9 13 ~8°, and 9 23 ~ 44°, 

where the new value for 6*12 lies within the 90% (and 95% and 99% and 99.73%) 
confidence level region of the MSW LMA solution [23] and the new value for 
#13 is still below the CHOOZ upper bound. Note that the new values for #23 is, 
in principle, the same as the old one. At the same time, the exact values of the 
charged lepton masses can be accommodated by choosing the values Y" a £ = 0.5, 
Y* = 1.8, and Yj = 1.0 for the order unity coefficients. Hence, our model is 
in prefect agreement with the MSW LMA solution and it can reproduce the 
realistic charged lepton mass spectrum. 



8 Implications for neutrinoless double /3-decay, astrophysics, and 
cosmology 



Assuming massive Majorana neutrinos, we analyze the implication for the 
prediction of the effective Majorana mass \{m)\ in neutrinoless double /3-decay 
(see, e.g., Refs. [75,76]) 



(m)\ = 



U la m »a 



a=l 



(73) 



99.73% C.L.: 0.24 < tan 2 12 < Oi 



26° < \0 12 \ < 43° [23]; best-fit: 



tan 2 #12 ~ 0.42 



1 012 1 ^ 33° [23] 



10 , 



C.L.: sin 2 26 2 z > 0.92 [9]; best-fit: sin 2 26> 23 ~ 1.00 



1 023 1 ^45.0° [8] 
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where the Z7i 's are first row matrix elements of the leptonic mixing matrix U 
and the m^s are the masses of the neutrino mass eigenstates. Inserting the 
expressions for the Ui a 's and the m^'s into Eq. (73), we obtain 



1 



(m) | = -m u [(V8 + e 4 cos 20 12 + e 2 ) cos 2 13 + 4e 4 sin 2 13 



+ 



V2 cos 29i2 cos 2 813 H — cos 2 6> 13 e 2 



cos 20 12 cos 2 0i3 + 2 sin 2 13 e 4 + C(e 8 



~ m u y/2 cos 2#i 2 cos 2 6*13 



(74) 



Note that the quantity | (m) | becomes equal to zero when On = f , i.e., when we 
have an exactly maximal solar mixing angle.|33 Using m u ~ 0.04 eV, On — 41°, 
and #13 ~ 4°, we find that | (m) \ ~ 0.007 eV, which is consistent with and below 
the phenomenological upper bound < 0.080 eV [78] for an inverted hier- 
archical neutrino mass spectrum. Our value for | (m) | is also below the exper- 
imental upper bound obtained from neutrinoless double /3-decay experiments 
( 76 Ge experiments) by the Heidelberg- Moscow collaboration [|(m}| < 0.35 eV 
(@ 90% C.L.)] [79] and by the IGEX collaboration [\(m)\ < (0.33 1.35) eV 
(@ 90% C.L.)] [80,81]. It would also be below the sensitivity of the planned 
neutrinoless double /3-decay experiments GENIUS, EXO, MAJORANA, and 
MOON, which is |(m)| ~ 0.01 eV. 



Furthermore, the sum of the neutrino masses 



(75) 



a=l 



where the m a 's are the physical masses of the neutrino mass eigenstates, is 
often used in astrophysics and cosmology. Inserting the expressions for the 
m a 's into Eq. (75), we obtain 



M = m u (V8 + e 2 + 2e 



m, 



2^+ 2 + 



128^ 



0(e 



10 \ 



2y/2- 



m,. 



(76) 



Again using m v ~ 0.04 eV, we find that M ~ 0.1 eV. There exist several upper 
bounds for this quantity from several different branches of astrophysics and 



11 A maximal solar mixing angle could, e.g., imply the presence of a superlight 
Dirac neutrino. For a treatment of naturally light Dirac neutrinos using the seesaw 
mechanism, see, e.g., Ref. [77]. 
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cosmology. For example, using the presently best available data from cosmic 
microwave background radiation and large scale structure measurements, we 
have the upper bound M < (2.5 -=- 3) eV (@ 95% C.L.) [82,83]. Other upper 
bounds of this quantity come, e.g., from cosmic microwave background radi- 
ation, galaxy clustering, and Lyman Alpha Forest measurements M < 4.2 eV 
(@ 95% C.L.) [84] and from the 2dF Galaxy Redshift Survey M < (1.8^2.2) eV 
(@ 95% C.L.) [85]. Thus, our value for the sum of neutrino masses is well be- 
low the upper bounds derived from astrophysics and cosmology. However, it 
should be possible to combine cosmic microwave background radiation data 
from the MAP/Planck satellite with Sloan Digital Sky Survey measurements, 
which could give an upper bound of M < 0.3 eV [86]. Then, our obtained 
value for the sum of the neutrino masses is rather close to this upper bound. 



9 Summary and conclusions 



In summary, we have presented a model based on flavor symmetries and 
higher- dimensional mass operators. This model is a modified and extended 
version of the model given in Ref. [49] and it naturally yields the mass matrix 
textures 



e 3 e 2 e 4 ^ 



e e 



e 3 e 2 



/ 



and 



-1 



1 t A 
-1 e A 



for the charged leptons and the neutrinos, respectively. Note that these tex- 
tures involve the same small expansion parameter e. Our old model [49] had 
a symmetric mass matrix for the charged leptons, whereas the new model 
has an asymmetric one. Thus, our new model predicts the realistic charged 
lepton mass spectrum (with the order unity coefficients Y* = 0.5, Y* = 1.8, 
and Yji = 1.0), an inverted hierarchical neutrino mass spectrum, a large (but 
not necessarily close to maximal) solar mixing angle #12, which is in excellent 
agreement with the MSW LMA solution of the solar neutrino problem, a small 
reactor mixing angle 13 , and an approximately maximal atmospheric mixing 
angle #23 (enforced by the flavor symmetries). Furthermore, these textures 
yield the MSW LMA solution and it follows from the model that 612 — \ — #13. 
The explicitly obtained values for the mixing angles (assuming no CV viola- 
tion, i.e., 5 = 0) are 



012 ~ 41°, 13 ~ 4°, and 9 23 ~ 44° 



for the ratio of the order unity coefficients Y^/Y^ = 1 and 



6 12 ~37°, 9 n ~8°, and 6 23 ~ 44°, 
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for the ratio of the order unity coefficients Y b e /Yj = 2. Both these sets of the 
leptonic mixing angles are in very good agreement with present experimen- 
tal data. In addition, our values are not in conflict with limits coming from 
neutrinoless double /3-decay, astrophysics, and cosmology. 

Finally, we have also (in the Appendix) presented a scheme for how to trans- 
form any given 3x3 unitary matrix to the "standard" parameterization form 
of the Particle Data Group. 
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A Transformation of any 3x3 unitary matrix to the "standard" 
parameterization form 



The "standard" parameterization form of a 3 x 3 unitary matrix according to 
the Particle Data Group [28] reads 



U = 



C2C3 S3C2 S2G 1<5 ^ 

-S3C1 — SiS2C^e^ C1C3 — S'liS^iSjse 1 ' 5 S1C2 
S1S3 — S^CiC^e 15 —S1C3 — S^^Cie 115 C1C2 j 



(A.l) 



where S a = sin# a , C a = cos9 a (for a = 1,2,3), and 5 is the physical CV 
violation phase. Here 9\ = #23, #2 = ^13, and 63 = 612 are the Euler (mixing) 
angles. Note that four of the entries of the matrix U in Eq. (A.l) are real, i.e., 
the 1-1, 1-2, 2-3, and 3-3 entries. Furthermore, the matrix U in Eq. (A.l) can 
be decomposed as follows: 
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U = 23 (9 23 )U 13 (9 13 ,5)0 12 ( 
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■S 13 e iS 


C 13 


) 


V 
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(A.2) 



where C a b = cos^, = sin^ (for a, 6 = 1, 2, 3), and a b(9 a b) is a rotation 
by an angle a & in the afe-plane. If 5 = 0, then L r i3(6'i3, 0) = 013(6*13). 

Here we will show how to transform any given 3x3 unitary matrix to the 
"standard" parameterization form. A 3 x 3 unitary matrix 



^11 U 12 U 13 ^ 



u = (U ab ) 



U 2 \ u 22 u 23 
U 3 i u 32 u 33 



(A.3) 



where a,b = 1, 2, 3, obeys WU = UW = 1 3 (9 unitarity conditions) and is 
therefore characterized by 9 real parameters, since a general 3x3 complex 
matrix is characterized by 18 real parameters (9 complex parameters). Out 
of these 9 parameters, 3 are Euler (mixing) angles and 6 are phase factors 
(phases) . Not all of the 6 phases enter into expressions for physically measur- 
able quantities. Only those phases are physical], which cannot be eliminated 
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In the transition probability formulas for neutrino oscillations (y a — > vp) 

_3 3 A _2 

P(y a -> up) = P a p = P a p(L) 



EE* 4 - 



a=l 6=1 



where a, (5 = e, //, r; L is the neutrino baseline length, E is the neutrino energy, 



J; 



ab 

a/3 



Ua a Up a U a bUp b are the amplitude parameters, and Am^ 



mr are 



the neutrino mass squared differences, we observe that the matrix elements of the 
leptonic mixing matrix U = (U aa ), where a = e, /i,t and a = 1,2, 3, only appear in 

jab 
'a/3- 

parameters 



the amplitude parameters J°~a. It is obvious from the definitions of the amplitude 



JaP = UaaUpaUabUpt, 

to see that they are invariant under phase transformations of the following kind 



Ka 



o il PaTJ 



where <p a ,ip a £ K are arbitrary parameters, i.e., J'^p = J®p. Thus, the transition 
probability formulas may only depend on phases in the matrix U, which cannot be 
absorbed by the above phase transformations (see Eq. (A. 4) for the corresponding 
matrix form). The number of such phases is equal to 1 (in the 3x3 case), i.e., the 
CP violation phase 5 [87,88]. 
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by the transformation 



where the matrices $^ = diag (e 1 ^, e 1 ^, e lipT ) and & u 
can always be put on the forms 



(A.4) 

diag(e ivi ,e i ^,e i ^ 3 ) 



(A.5) 
(A.6) 



such that det$£/ = 1 and det<Ev = 1, where <Ey = diag (e i¥, <=' , e il/ V , e'^' ) and 
<JV = diag (e i</v , e'^ 2 ', e 1 ^'). This means that cp a = (fe + </V (« = e, /i, r) and 
V^a = + </V (a = 1, 2, 3), where </v + </V + yv = and </?i/ + </? 2 ' + V?3' = 0. 
Thus, we have 



U^U' = e^-^QeU®^ = e i(f <f> e U<S> v/ , 
where ip = ^ — y^, and on matrix form, we find that 

Vne^'-^') C/ise 1 ^'-^') C/^e^'-^')^ 
C/ 21 e i( ^'-^i') [/^eHw-^') [/^eK^'-vsO 
[/gjeKiPr'-viO [^ae'^r'-w) C^e 1 ^'-^') 



(A.7) 



C/' = e 



V 1 



(A.8) 



Assuming now that the matrix U' in Eq. (A.8) is on the "standard" param- 
eterization form as displayed in Eq. (A.l), we identify from the 1-1, 1-2, 1-3, 
2-3, and 3-3 entries that 
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= e iv U u e 
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-<Pl>) — 




e i{ 
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= e iip U 12 e 
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e i{ 




= e^U 13 e 


(W 


-<P3>) — 


\u 13 


e i( 


S\C 2 
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argf/ 2 3+y+¥> M /-¥v) 
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(A.9) 
(A.10) 
(A.11) 
(A.12) 
(A.13) 



Taking the imaginary parts of Eqs. (A.9), (A. 10), (A.12), and (A.13), we arrive 
at 



= arg Uu + (p + (p e , - ip v , (A.14) 

= argL r i2 + (p + (fie' - <f2', (A.15) 

= arg£/ 23 + if + </y - ipy, (A. 16) 

= arg[/ 33 + <p + <p T > - </? 3 /, (A. 17) 

which together with the two conditions 
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W + W + W = 0, 

(fill + If 2' + f 3 > = 



(A.18) 
(A.19) 



make up a system of equations. Note that the overall phase f can be absorbed 
into the argf/ ab 's, i.e., &rgU a b + f — > arg£/ a 5. The system of equations includes 
six equations and six unknown quantities, i.e., f a > (a = e,/j,,r) and <f a > (a = 
1, 2, 3), and thus, it has a unique solution. This solution is 



fe 



f2' 



(-2argE/n - 2argE/i 2 - arg£/ 23 - arg£/ 33 ) , 



W = - (argC/n + arg£/ 12 - arg£/ 23 + 2arg£/ 33 ) , 
W = -yv - </V> 

¥>i' = ^ (argC/n - 2argC/i 2 - arg£/ 23 - arg£/ 33 ) , 



'-2 arg U u + arg C/i 2 - arg U 23 - arg C/ 33 ) , 



<f 3' = -fV - f>2'- 

Inserting the solution into Eqs. (A. 9) - (A. 13), we obtain 



(A.20) 

(A.21) 
(A.22) 
(A.23) 

(A.24) 
(A.25) 



C2C3 — 
S3C2 = 
S 2 e~ i5 = 
SiC 2 = 
C\C 2 = 



Uu\, 
Uu\, 

TJ |„i(argC/i3-argC/n-argl/i2-argl/23~arg(733) 

U23I 
U 33 \. 



(A.26) 
(A.27) 
(A.28) 
(A.29) 
(A.30) 



Thus, it follows from Eq. (A.28) by taking the real and imaginary parts, 
respectively, that 



S2 = |^i3| =>• #13 = &2 — arcsinS^ = arcsin \U\ 3 \, 
-5 = - arg U u - arg U 12 + arg U 13 - arg U 23 - arg U 33 . 

Dividing Eq. (A.27) by Eq. (A.26), we find that 



5s 
C3 



U 



12 



u 



11 



? 12 = 9 3 = arctan ^ = arctan 
^3 



U 



12 



U- 



11 



and similarly, dividing Eq. (A.29) by Eq. (A.30), we find that 



u. 



23 



17: 



33 



9 23 = 9 1 = arctan = arctan 



U- 



23 



33 



(A.31) 
(A.32) 



(A.33) 



(A.34) 



In summary, the parameters of the "standard" parameterization form is given 

as 



36 



e 



12 



arctan 



U 



12 



6 'i3 = arcsin | C/x3 1 , 

C/23 







23 



arctan 



t/33 

(5 = argL r n + argC/12 



arg Z7i3 + arg U23 + arg C/ ; 



(A.35) 
(A.36) 
(A.37) 
(A.38) 



in terms of the entries of any 3x3 unitary matrix U = (U a b)- Note that 
these parameters are uniquely determined by the entries in the first row and 
third column (Un, U\ 3 , C/ 23 , and U 33 ) only as well as the overall phase 
if. The other entries are completely restricted by and follow directly from the 
unitarity conditions (WU = UW = I3). 
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